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Abstract —Quaternion has many excellent features to 
describe rotational motion of rigid body, so it can be used to 
describe the process of tool axis vector. By making unit tool axis 
attitude points correspond with the quaternion spatial points, 
the optimization of two axis vector can be converted into the 
optimization of the quaternion spots in quaternion space. In 
quaternion space, through the use of B-spline curve and global 
interpolation method, the quaternion spots can be interpolated 
and optimized. In this way, it is able to achieve that fit tool axis 
vector, optimize machining process of tool axis trajectory, and 
make the cutter path smoother. 

Index Terms — Quaternion, B-spline curve, global 
interpolation method, cutter axis vector, smoother 


I. INTRODUCTION 

Be different from three axis NC machining, cutter axis 
vector can be changed flexibly in five axis machining, this 
difference is caused by the two axes of rotation. Because of 
the two axis of rotation, the five axis CNC machine can 
control tools to finish works at arbitrary attitude, especially 
in the complex curved surface machining at high speed and 
precision, have demonstrated irreplaceable advantages. 
However, if two rotating shaft became more flexibility, it 
also increase the complexity of the structure of the machine 
tool, make it difficult to visualize the actual processing of 
the machine tool movement , and increase the difficulty to 
control the tool posture vector. In addition, same with the 
feeding speed, abilities of tool axis rotation also have 
limitation. 


Literature [2] was first used C-space to search the 
interference free tool axis vector of each cutter location, and 
in each cutter axis vector C-space range, calculate the global 
optimal cutter axis vector according to the principle of 
minimum rotation of adjacent cutter. But the problem is that, 
on the one hand, the method have too much iterative, result 
in serious decline of efficiency, and on the other hand, the 
method still cannot guarantee a smooth transition between 
adjacent cutter axis vectors. The literature on [3] and [4] 
respectively used the method of quaternion linear 
interpolation to smooth cutter axis vector between two key 
cutter shafts. All these methods are just applicable to a 
smooth transition between the two key tool axis vectors, and 
cannot handle the transition between multiple key tool axis 
vectors. The literature on [5] research literature only one 
cutter axis vector, without considering the problem of 
smooth transition of global tool axis vector. 

In the actual machining, CNC system always based on 
CAM system to generate the tool posture vector, and the tool 
vector points can be regarded as unit vector points. Fig.l 
shows the process of tool axis vector optimization. 



The sudden changes of the tool axis vector are easy to cause 
the vibration of the machine, result in the declining in the 
quality of processing. Therefore, in five axis NC machining 
process, the tool axis vector should be avoided having 
sudden changes, which include its rotation direction and 
rotation speed. 

In process of manufacturing, feed speed and rotational 
speed are fitted on each other. When tool axis vector rotation 
is large between two adjacent cutter contacts, in order to 
meet the requirement of the rotating speed, the feed speed 
should be reduced to relatively lower. It likes that the cutter 
is staying in one position for a long time, and leaving cut pit 
more obvious in the part surface, affecting the quality of 
parts. So, it is necessary to smooth the tool axis vector. 
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Fig.l Flowchart of the tool axis vector optimization based on 
quaternion 

II. QUATERNING BACKGROUND 

The so-called quaternion is defined by a number of unit 1 
and three imaginary unit i, j and k. The general form 

is A = (A 0 ,A l ,A 2 ,A 3 ) = Aq + \i + \j + ^k , and the 

unit 1, i, j, k can be regarded as the unit total in 
four-dimensional space vector. Any of the quaternion 
numbers can be represented by points or vectors in space. 
Quaternion could be well used to express rotation, and 
provides a very convenient way for the study of rigid body 
motion. Quaternion operation have many properties, and the 
properties of the operations were used in this article are:® 
When a quaternion number is multiplied by a scalar, all 
elements should be multiplied by this scalar, in particular, 
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the negative of quaternion number is 
— A = — Aq — A^i — A^j — A^k .(2) The quaternion element 
of product is still quaternion number, the algorithm is that 
i © i = —1 , j o j = -1 , k o k = -1 , i o j = -j o / = £ , 

koi = -iok = j,jok = -koj = i. 

III. QUATERNION OPTIMI Z ATION TOOL POSTURE 

Quaternion provides a very convenient way for the study of 
rigid body motion, because of the quaternion, the rigid body 
movements can be described conveniently and intuitively. 
For a series of tool posture vector are generated by CAM 
such as p v p m ,and they can be expressed by 

quaternion as p i =Q i p l Q7 l ( i =1,2 ’ ... ’ 

m ) . (1) 

Here, Q.=A V + A^i + A^j + A^k , p x =i, and Q i =E. 

Given p t = p a i + p i2 j + p i3 k , and by compared with (1), 
the equations can be obtained as: 

Af) + \ ~ ^2 ~~ ^3 = Pi\ 

< 2(A v A 3 + \A, 2 ) = p i2 .(2) 

2(A l A 3 +A 0 A 2 ) = p i3 

By the formula (2), the expression of quaternion points 
<2- will be expressed as follows: 

Q t = ±-^(l + co i )(cos20 i +sin2^.)x(/+ ^ i2 j+ ^ i3 k) 
V2 1 + Pa 1 + Pn 

.(3) 

Here 9 i are angles around the axis of rotation. Since 
formula (3) it can obtain a series of quaternion points: Q x , 
Q 2 ,Q 3 ,.•••, Q m • As shown in Fig. 2, the rotation of tool 
posture points is described in quaternion space as follows: 



Fig. 2 The tool posture points rotation in quaternion space 

IV. INTERPOLATE THE QUATERNION POINT BY 
B-SPLINE 

Now it is going to use the global interpolation method to 
interpolate quaternion points Q x , Q 2 , Q 3 ,...., Q m and get the 


I a-a , I 

First, it can get V- from V- = V M H-, where 


d — ^1 Q. — Q i X I , i= ( 2,3 ’ ... ’ m-1) , then according 


to: 


i=i 

M 0 =... = U p =0,U t _ p =... = U, =1 

l J+P~ l 

U j+P = ~ 'Zv i ,j = l,2,...,m-p 

P i=j 


the node of 


vector U={ u 0 ,u x , ., U t }will be obtained. According to 

U, 


B UP (y,) = - 


0, otherwise 


B,, P (v,) = • B. p _, (v,) + V ' • B lw (v,.) 


w. —u 

l+p l 


U i+p +1 *L+1 


the matrix equations can be expressed as : B j (y.) • q { = • 
Then, control points q t can be obtained and B-spline curves 


in 

of quaternion be expressed as: <2/ = ^ „(v,)g, • As 

i=i 

shown in Fig.5, the quaternion points <2- are fitted by 
B-spline curve: 



Fig.3 B-spline fitting quaternion points 


V. FITTING TOOL POSTURE VECTOR 

When get the B spline curves of quaternion, it can 
interpolate and optimized the discrete tool posture vector in 
quaternion space, so as to get the tool motion curve smoother. 
Now it starts to optimize tool orientation as following: 

m 

When the quaternion point as Q j = ^ B j p (v ; )q t was 

i=i 

obtained, and according to the expression p t = Q i p 1 Q i 1 , 
the other expression of tool posture can be derived as follows: 

m m 

p t == 

i=l J =l 


quaternion curve Q i 
control points. 


m 

5 p ( V i )q i , here q { are unknown 

i=i 


zm ✓-» ✓-» 

! ' +>1 c i+/ 
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p p p p 

T Tj/- v* Ci Cj c i Cj i 

Let w = 2 ^ 2^~ q$fl\ - ip q t qfa > and too1 

i+J Ci +j k Ck 

orientation curve expression could be got 

2m 

as.p, = 2X, | ,(v,).W'. 

k =1 

Based on the algorithm stated above, the tool orientation 
curve will be got and described in the unit sphere as 
illustrated in the following Fig.4, we can see that tool 
orientation curve fitting is performed based on quaternion 
points. 



Fig.4 The tool orientation curve in the unit sphere 


VI. ANALYSIS OF THE ALGORITHM 

In order to test the validity of this algorithm, the following 
work is on a track of the actual simulation. Fig.5 compared 
the normal vector of tool posture points between before and 
after optimized as following: 



(a), the normal vector before optimization 



smoother after optimized. This confirms that the tool axis 
vector optimization method in this article is effective. 

VII. Conclusion 

For five-axis machining, if the angles of cutter shaft 
between two adjacent cutter locations are too large, or the 
cutter axis vectors are uneven, all these would make the 
cutter axis influence the surface quality. In this article, 
through the use of quaternion on the rigid body rotation and 
convenient description, the tool posture can be converted 
into quaternion and transform the discrete tool posture 
interpolation problem into the curve design interpolation 
problem. In this way, the tool vector can be optimized and 
the quality of machining can be improved. 
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(b). the normal vector after optimization 

Fig.5 The normal vector of the tool posture point 

Fig.5 (a) shows the normal vector before optimization, and 
Fig.5 (b) is the normal vector after optimization. It is evident 
from these figures that the cutter axis vector becomes 
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